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We simulate Kn scattering in s wave and p wave for both isospins I — 1/2, 3/2 using quark- 
antiquark and meson-meson interpolating fields. We extract the elastic phase shifts S at several 
values of the Ktt relative momenta. The resulting phases exhibit qualitative agreement with the 
experimental phases in all four channels. We express the s wave phase shifts near threshold in terms 
of the scattering length and the effective range. Our Kn system has zero total momentum and is 
simulated on a single ensemble with two dynamical quarks, so results apply for ~ 266 MeV 
and mpc — 552 MeV in our simulation. The backtracking contractions in both 7=1/2 channels 
are handled by the use of Laplacian-Heavyside smeared quarks within the distillation method. 
Elastic phases are extracted from the energy levels using Liischer's relations. In all four channels 
we observe the expected K(n)ir(— n) scattering states, which are shifted due to the interaction. In 
both attractive 7=1/2 channels we observe additional states that are related to resonances; we 
attribute them to if,* (1430) in s wave and A"* (892), if* (1410) and if* (1680) in p wave. 

PACS numbers: 11. 15. Ha, 12.38. Gc 



I. INTRODUCTION 



Lattice QCD (LQCD) provides an approach for ab ini- 
tio calculations of hadron properties. Only in recent 
years, the tools are becoming efficient enough to study 
strong interactions between hadrons and strong decays 
of hadronic resonances. Due to the finite spatial volume, 
the spectral density of scattering processes is intrinsically 
discrete, and thus one has to infer scattering information 
from the observed energy levels. In current simulations 
these levels are scarce, and with available methods, phase 
shifts S can only be determined at a limited number of 
values of the invariant mass y/s. Given these limitations, 
it is promising that recently several lattice studies [3-0| 
determined the tttt scattering phase shift in p wave with 
7 = 1, where the p meson dominates. The related 7T7t 
scattering in s wave with 7 = 2, which does not require 
backtracking contractions, has been thoroughly explored 
on the lattice recently |ci-[l(j|. 

Continuing along that path, we study here the 7f7r 
system in s wave and p wave, for both, 7 = 1/2 and the 
exotic 7 = 3/2 channels. Whereas the I — 1/2 p wave 
is dominated by the well established and narrow reso- 
nance if*, the experimental 7 = 1/2 s-wave phase shows 
a very broad structure without clear resonance signal be- 
low 1 GeV. The discussion whether this should be inter- 
preted as a wide resonance 7f q (800) (also called k) with a 
width almost as large as its mass is continuing. It would 
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fit into the + multiplct of scalar mesons (together with 
the partner states /o, also called a, and oq). 



A. Experiments 



The experiments provide the magnitude and phase of 
the scattering amplitude 1 Tj = |T/|e 1( ^ as a function of 
the 7f7r invariant mass \fs. In the elastic region, both 
|T/| and <$>\ = 5{ are related to the elastic scattering 
phase Si 
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In early experiments, 7f7r scattering amplitudes were 
derived from analysis of Kp — > Kirn and Kp — ¥ Ktt A 
[llT - flfjj . Even today, the most accurate data on scat- 
tering amplitudes is based on Estabrooks ll| and As- 
ton [12j . and we use these two for comparison with our 
lattice results. Estabrooks ll| measured both processes 
and was the only experiment that was able to disentangle 
7=1/2 and 7 = 3/2 channels. Their phases 5\~^ 2,3 ^ 2 
are plotted in Fig. [21 and all the presented points from 
Ref. [ll| are in the elastic region. 

Aston [l2| considered only 7fp — > Kirn, so they were 
able to provide only the sum Ti = T^ 2 + \T^ 2 = 
\Ti\^ 1 . The phases ^?J^ X are believed to be small (which 
was explicitly confirmed by Estabrooks), so we compare 
to their dj~^ 2 ~ 4>i=i in Fig. [2j We also compare to 



Our T{ is denoted by a( / \J2t + 1 in experiments llll . [l2Tl . 
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their S^^ 2 , which is obtained taking their To [l2[ and 

subtracting 2 T 7=3/2 [nj]. The (blue) stars in Fig. [2] cor- 
respond to Aston's phases which are fully in the elastic 
region, i.e., the amplitudes lie on the unitary Argand 
circle with radius |2T/ — i| = 1. The (green) crosses 
present measured phases b\ = <3?f from the same exper- 
iment for which the amplitudes are almost elastic, i.e., 
the radius of the Argand circle is allowed to be in the 
range 0.85 < |2T/-i| < 1.15. 

The kinematics of K 7r-scattering, the analytical struc- 
ture of the partial wave amplitudes and the experimental 
results until the late 1970s have been reviewed in Rcf. 

Newer results on Ktt scattering were derived from D — > 



KtTTT 



- 4201 ] and B — > D*K and sequential D* decay 
[2l| - |24| . but none of these studies performs the isospin 
decomposition of the amplitudes. 



B. Theory: Continuum 

In particular the existence and parameters of the scalar 
if q (800) and its partner scalar nonet states have been a 
continuing source of discussion. 

Continuum calculations have been based on unitarized 
quark models [25T - [29j , chiral perturbation theory (ChPT) 
llHH and unitarized ChPT [HH. 

Unitarized ChPT expansions have also been used to 
study finite volume effects [40l - [43j in order to understand 
what features to expect from lattice calculations, partic- 
ularly for the scalar channel. 

The scattering amplitudes were parametrized in the 
most general way allowed by quantum field theory in 
Refs. 0, |45| according to the Roy-Steiner approach; 
the position of the if * (800) and if* (892) poles in the 
complex plane were then derived by using the experi- 
mental knowledge of Ktt scattering amplitudes at high 

Examples of further analytical studies related to the 
Kq (800) are given in Refs. |46l448| and references therein. 



considered p wave Ktt phase shifts, and we aim to deter- 
mine them here. 

The extraction of the s-wave phase shifts near the 
threshold was mainly focused on the 1 — 3/2 chan- 
nel K + tt + , since it does not require the evaluation of 
challenging backtracking contractions. The scattering 
length was determined from the finite volume energy 
shift AE = E — m n — rriK utilizing Liischer's formula 
[49l | . The quenched simulations [50|, HH were followed by 
NPLQCD using 2+1 staggered sea quarks and domain- 
wall valence quarks [13, by PACS-CS using 2+1 Wil- 
son sea and valence quarks on V = 32 3 x 64 [53| and 
by Fu using 2+1 staggered sea and valence quarks [lij]. 
We compare all results as a function of m T in Fig. [5] 
by showing the ratio o-o/^Ktt with the reduced mass 
Ma'tt = m-^mK / (Tn n + mx), since the quantity an/ ' IXk-w 
is i ndep endent of m,^ in lowest-order ChPT. ChPT 
[53l |54| or mixed ChPT [52[ is used to extrapolate the 
results derived at higher pion masses down to the physi- 
cal point. The ChPT expansions for the Ktt system are 
considered in Refs. [5514601 ]. 

The I = 1/2, I = channel involves also challenging 

1—1/2 

backtracking contractions. The scattering length a 
was determined in the quenched simulation 15 1||. then in 
the dynamical studies by PACS-CS [HH and Fu |54| , men- 
tioned already above. The results are compiled in Fig [5] 
NPLQCD [H^l extracted Oq -1 ^ 2 in the chiral limit only 
indirectly through the knowledge of low energy constants, 
without actually simulating the 1 = 1/2 contractions. 

1 /2 

The extraction of the phase shift 8 from the first 
excited energy state was actually done in a simulation 
with dynamical staggered fermions in Ref. jolj . How- 
ever, note that the ground state in such a simulation 
corresponds to the staggered taste if 57^5 while differ- 
ent unphysical tastes K^tt^ with 6^5 [62| - [64j 3 are ex- 
pected to contribute to excited states. Therefore, the 
phase shifts extracted from the excited states in such a 
simulation may correspond merely to staggered artifacts 
rather than physics of Ktt scattering. 



C. Theory: Previous lattice studies 



Up to now, lattice simulations of Ktt scattering have 

close to 



extracted only the s-wave phase shifts 5 



1=3/2,1/2 
1=0 



the threshold, which are commonly expressed in terms 

j 1/2 3/2 

of the scattering lengths a ~ , defined in Eq. CT~ 



Lattice simulations have not yet extracted the Ktt s-wave 
phase shifts away from threshold and that is one of the 
purposes of the present work. Simulations have also not 



The indirect lattice determination of the Ktt s-wave 
phase shifts in the I = 1/2 channel was addressed 
through the simulations of the scalar semilcptonic K — > tt 
form factor fo in Ref. (6|| . 

To summarize, there has been no direct simulation of 
the s wave if7r phase shifts away from threshold, and 
no simulation of the p wave Ktt phase shifts. These are 
addressed in the present work. 



2 Instead of subtracting the measured Tq ^ 2 [Tll | we subtract the 
effective-range formula fit through this data. 



3 The Ktt interpolator can be projected to desired taste of both 
mesons, while the su interpolators used in fsitl inevitably couple 
to all tastes as shown in Refs. [62|464| | . so the variational analysis 
is expected to render K^tti, as excited states. 
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II. ANALYSIS TOOLS 
A. Energy levels and phase shift 

In LQCD one determines Euclidean correlation func- 
tions, in the simplest case those of products of two inter- 
polators Oij with the quantum numbers of the hadronic 
channel at some Euclidean time distance. For finite vol- 
umes the spectral function is no longer continuous and 
thus all information has to be derived from the discrete 
energy levels. These are hidden in the spectral represen- 
tation of the correlation function for a set of interpolators 
0(t) with the same quantum numbers 



C jk (t) = (O i (t)O fe (0)t) = yVo^OjWe- 



E„t, 



\O k (0) 



t\ 



(2) 

The state-of-the-art method to recover the low lying en- 
ergy levels E n is the variational method [66l-l69j . The gen- 
eralized eigenvalue problem C(t)u n (t) = X n (t)C(to)u n {t) 
then disentangles the eigenstates n making it possible to 
obtain energy levels from the exponential decay of the 
eigenvalues 



while the effective energies 



E n (t) = log 



An(*+1) 



(3) 



(4) 



render E n (t) ~ E n at large t. In actual calculations it 
is not possible to have a complete set of interpolators al- 
lowing one to represent the physical eigenstates. One is 
limited to a reasonable subset. Also, the statistical qual- 
ity of Cjkif) is an issue. The reliability of the obtained 
energy levels decreases for higher \n), the ground state 
being the most reliable one. 

Effective mass plots E n (t) [see Eq. Q] are used only 
to estimate the fit range for the exponential fits to the 
eigenvalues. The energy values are extracted using cor- 
related fits of \ n (t) to one or two exponentials. When 
using two-exponential fits starting at small t, we verify 
that the extracted levels agree with results obtained from 
one-exponential fits starting at larger t. 

The information on the Kir scattering is contained in 
the scattering amplitudes T/(s) in Eq. ([T]). Here we 
concentrate on projections to isospin / = 1/2,3/2 and 
partial waves I = 0, 1. We choose the total 3-momentum 
P of the Kir system to be zero, so the lattice frame repre- 
sents also the center-of-momentum frame in our current 
simulation: 

P = Pn+PK = , P7r = -pa- = p* , s = E 2 -P 2 = E 2 

(5) 

We measure the energy E of the interacting Kn system 
in finite volume 



E = = \(p 



P 



w\2 



= Vp tt 



*2 



K ' 

(6) 



(here p K and p" denote 4-momenta) which allows us to 
extract the momentum p* = |p*| and the dimensionless 
q via 



V = 



(m K + raO 2 ] [s - (m K - m^f 
4s 



q = P 27 



(7) 

Assuming that the strong interaction of Kit is localized 
to r < R, the extracted p* represents the momenta of it 
and K in the outer region r > R. The resulting momenta 
p* or q can be related to the scattering phase shift in the 
elastic region 0, E3, [lO, [ill 



tan 6(q) 



7T 3 / 2 C 



2 00 (i;<7 2 ) 



for P = 



(8) 



This is Liischer's relation [70( between the energy levels 
in finite volume [which enter via Eq. (J7|)] and the phase 
shifts. The generalized zeta function Zi m is given in Ref. 

M- 

In the discussion of our results, we also apply the com- 
bination 



pl(s)=W2— cot Sf(p*) 



(9) 



which vanishes at the position of the resonance 8(sr) = 
t;. The variable p\(s) provides a convenient parametriza- 
tion of the elastic partial wave near threshold, 



r spl( S ) = \ + kr I eP * 2 + 0(p* 4 ) 



(10) 



where a is the scattering length and r the effective range. 
Near a relativistic Breit-Wigner resonance in the elastic 
region we may write the partial wave amplitude [72| 



T 



-Vir(s) 



s - s R + i v / ir(s) 



e lS '^ s) sin6j(s) . (11) 



where sr = m 2 R denotes the resonance position and T is 
the decay width. Considering the threshold behavior we 
can define 



and get 



r(s ) = [ PJ — 7 



/ir(s) cot 6(q) = pt(s) 7 = s R 

Pei s ) = ~( S R - s ) ■ 

7 



(12) 



(13) 



For a resonance in the I = 1 channel one has 7 = g 2 / (6tt), 
defining the coupling constant g (e.g., gK*Kn)- 

Naively, one would expect that in simulations with dy- 
namical quarks, all possible intermediate states should 
contribute. In the correlation function of a p meson with 
a quark-antiquark interpolator, one thus would expect to 
find signals of the 7T7t intermediate state (in p wave). It 
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turned out that this appears not to be the case, most 
likely due to weak coupling of the corresponding lattice 
interpolators. One had to include tttt interpolators ex- 
plicitly in order to obtain energy levels representing scat- 
tering This motivates us to include Kir and some 
other meson-meson interpolators in the present study. 

In the study of the nn scattering pj, we also used in- 
terpolators with total momentum P ^ 0, which allowed 
us to find phase shifts at more values of s = E 2 — P 2 . 
The corresponding phase shift relation (JSJ for scattering 
of particles of equal mass mi = mi and total momentum 
P 7^ has been derived in Refs. [z^, [lH- For systems 
with two mesons of different mass mi ^ mi and P ^ 0, 
the formalism has been extended in Refs. (42l. l75T - [77j . 

In our present study, we consider only the case P = 
Ptt + V>k = 0, and so the original phase shift relation (0) 
applies [201. We consider the irreducible representa- 
tion of Oh to extract the s wave and the T-[~ irreducible 
representation to extract the p wave. A± will also con- 
tain the admixture of I = 4 and higher partial waves, and 
TT will contain admixture of I = 3 (and higher waves) 
[70j . but those are expected to be small, and we neglect 
the effect of such mixing. In contrast to that, the sim- 
ulations of Ktt scattering with P / would involve an 
additional complication since the partial waves for even 
and odd £ can mix in the same irreducible representa- 
tion [13, l75l - [77j . We avoid this additional complication 
by taking P = in the present work. 

It is important to provide a large enough set of in- 
terpolators for a good representation of the lowest phys- 
ical states in the generalized eigenvalue analysis. The 
interpolators used for the four channels studied (isospin 
1/2 and 3/2, s wave and p wave) are given in Appendix 
IA1 In addition to several qq operators 0^ q we also in- 
clude several meson-meson operators Q MM (Ktt, K* p, 
K\ ai ), in total up to eight for the s wave and up to six 
for the p wave. We include for example ir(Q)K(Q) and 

J2i=x,y, z i 1T (Pi) K (~Pi) + n(-Pi) K (Pi)] for s wave ; and 

n(p z )K(-p z ) - ir(-p z )K(p z ) for p wave (p, = ^e<). 
Here, each meson is separately projected to definite mo- 
mentum, given in parentheses in units of 2tt/L. The Wick 
contractions are provided in Appendix [B] 



B. Lattice simulation 



NlxNr 13 a[fm] L[fm] #configs m^MeV] m K [MeV ] 
16 3 x 32 7.1 0.1239(13) 1.98 280 266(3)(3) 552(2) (6) 



TABLE I: Configurations used for the current study have two 
light dynamical flavors. Nl and Nt denote the number of 
lattice points in spatial and time directions. For more details 
see Q. 

parameters used for the simulation along with the num- 
ber of (approximately independent) gauge configurations 
used, the lattice spacing, volume and the pseudoscalar 
masses m,! , which are the most relevant for Kir scat- 
tering (for details, see Ref. [l[). 

The s quark is included only as a valence quark in 
the hadron propagators. To determine the strange quark 
hopping parameter k s , we calculated the connected part 
of the <f> meson which is expected to be almost exclu- 
sively ss. The tuning was done on sources in a single 
time slice and we obtained n s = 0.12610. Using our 
complete set of perambulators calculated with this value 
of k s we again determine the mass of the <p meson on 
the full data set and obtain m 1 ^ = 1015.8 ± 10.8 MeV 
which has to be compared to the experimental mass 
m e * p = 1019.455 ± 0.020 MeV. This k s corresponds to 
the kaon mass provided in Table H] 

The sea and valence quarks obey periodic boundary 
conditions in space. The gauge field obeys periodic 
boundary conditions, while the sea quarks satisfy an- 
tiperiodic boundary conditions in time. We compute and 
combine valence quark propagators with both antiperi- 
odic and periodic boundary conditions. This effectively 
extends the time direction to 2Nt = 64 by combining 
the periodic propagator Mp 1 and antipcriodic propaga- 
tor M^ 1 (see for example [10,[ll[). All results in this pa- 
per have been obtained using the so-called "P+A" prop- 
agators 

M -i (t t ^ = [\[Mp 1 {t f ,t i ) + Mj l {t f ,t i )\ t f >U, 
p +A V f, i) \i[ M p\tj,U) - MxHtfM t, < t, . 

(14) 

For further discussion of this method, see Ref. [l[ . 



Like in Ref. [l| we use configurations generated for 
the study of reweighting techniques [78|, [79[ kindly pro- 
vided by the authors. The action used to generate the 
gauge configurations containing rif = 2 flavors of mass- 
degenerate light quarks is a tree level improved Wilson- 
Clover action with gauge links smeared using one level of 
normalized hypercubic smearing (nHYP smearing). The 
practical advantage of rif = 2 flavor simulation for Ktt 
scattering (with respect to n/ = 2 + 1) is that there is no 
Krj scattering state, so the inelastic threshold is higher 
(discussion in Sec. IIII B 1[) . The valence u/d quarks have 
the same mass as the sea u/d quarks. Table U lists the 



C. Propagators, contractions and distillation 
method 

The hadron correlation functions are constructed by 
combining quark propagators derived on the gauge field 
configurations. For the meson-meson interpolators used 
here, the Wick contractions lead also to backtracking con- 
tributions, for example depicted in the box diagram of 
Fig. [7J:. For a statistically reliable inclusion one thus 
needs all-to-all methods. The distillation method pro- 
posed in Ref. [82| provides these capabilities. It is based 
on separable quark smearing operators, i.e., a truncated 
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spectral representation of the unit operator in terms of 
the eigenvectors of the spatial lattice Laplacian. The 
quarks q s in the simulation are smeared according to dif- 
ferent smearing widths s [lj 

N v 

q s =^V fe V fe )t q (15) 
fc=l 

N v = 96, 64, or 32 
for s = n (narrow) , m (middle) , or w (wide) , 

rendered by the different number N v of lowest eigenvec- 
tors u' fe ' incorporated in the sum. The technique effec- 
tively replaces the quark propagators G(x — > y) by prop- 
agation from one source to another, so-called perambu- 
lators r(i —¥ j). This allows high flexibility in the Dirac 
and momentum structure of the hadron interpolators. 

We proceed as follows. First, the gauge links are four- 
dimensional normalized hypercubic smeared [83[ with the 
same parameters used for generating the gauge configu- 
rations: (ai, 02,0:3) = (0.75,0.6,0.3). On each gauge 
configuration, we calculate the lowest 96 eigenvectors of 
the lattice Laplacian on every time slice. For the calcu- 
lation of the eigenmodes, we use the PRIMME package 
[Hj|. For the determination of the quark propagators 
we use the df l_sap ecr algorithm provided in Liischer's 
DD-HMC package [83, |86[- Due to the large number 
of sources necessary for the distillation approach, an in- 
verter employing low-mode deflation techniques is espe- 
cially well-suited. 

Statistical errors are determined with a single elimina- 
tion jackknife procedure throughout. When extracting 
energy levels, we properly account for correlation in Eu- 
clidean time t by estimating the full covariance matrix in 
the given fit interval. For the covariance matrix, we use 
a jackknife estimate which is calculated on the ensemble 
average only. 



III. RESULTS 
A. Pion, kaon and dispersion relations 

Since we consider scattering of tt and K, we need 
their masses and their separate (noninteracting) energies 
E-Kji{p). For 7r and K (J p = 0~) we use the six inter- 
polators given in Eq. (|A1[) . with three smearing widths 
for each of the two Dirac structures. Their masses and 
energies are extracted from the variational analysis of the 
6x6 correlation matrix and listed in Table ILT1 

We find that the lattice energies i? Wi ^-(p=^n) agree 
for n 2 < 2 with the continuum dispersion relation within 
the error (see Table ILTl 4 ). So we use the dispersion re- 



lation E V) k(p*) = y m nK+P* 2 w ith m -K.K fixed to 
am T = 0.1673 and arriK = 0.3466 throughout the anal- 
ysis. This means that p* 2 and q 2 are extracted from the 
lattice energy E using Eqs. (|6I7|) . and the resulting q 2 is 
used to get the phase shift from Eq. (jHJ). 

As a cross-check, we have verified that our resulting 
phase shifts obtained in this way agree with the phase 
shifts obtained from the energy shifts AE n = E n — E n — 
Ek of interacting Kir with respect to the nearest non- 
interacting scattering level. These values are determined 
using the ratio X n (t)/[K(^) (*) ^K(n K )(t)} oc e~ A - E ™ t of 
the eigenvalues A„ (t) for interacting Ktt and the eigenval- 
ues for noninteracting if(njf) and 7r(n ff ) with momenta 

2ir n 

~ n K,TT- 

B. Kir scattering 

Before addressing the details of the analysis and the 
results for separate channels below, let us compare the 
main features of all four channels (s wave and p wave 
in I = 1/2, 3/2). The resulting energy levels for the 
Kir system arc presented in terms of the effective energy 
E{t)a [see Eq. ((4])] in Fig. [1] These levels correspond to 
the preferred interpolator choices listed in Table [TTTl The 
horizontal broken lines show the energies E = Ek + E^ 
of the noninteracting states K(n)Tr(—n) as measured on 
our lattice with p* = ^ -Jn. 

The resulting spectrum agrees with the expectations 
for the respective channels: there is a scattering state 
A'(0)7r(0) in ,s wave (black circles), which is below txik + 
m-n in the attractive 7 = 1/2 channel and above mK + m n 
in the repulsive 7 = 3/2 channel. There is no scattering 
state A(0)7r(0) in p wave due to nonvanishing orbital 
momentum. The scattering state A(l)7r(— 1) (green cir- 
cles at Ea ~ 0.95) is observed in all four channels, and 
its signal is nicer for I = 3/2 (with / = 0,1), since 
channels with maximal isospin I = 3/2 do not involve 
backtracking contractions. In addition to the scattering 
states, which lie close to noninteracting levels given by 
the dashed lines, there are additional states in the attrac- 
tive 1=1/2 channels (red and pink levels). While these 
states are of course also shifted with regard to the reso- 
nance position, we refer to these additional states as "re- 
lated to" the respective resonance. In this language, the 
additional state in the s wave is probably related to the 
scalar resonance Aq(1430), while the additional states in 
p wave are expected to be related to vector resonances 
A*(892), A*(1410) and A*(1680), respectively. 

The results for phase shifts S^(s), which are based on 
the energy levels of Fig. [I] are presented in Table IIIII 
They are compared to the experimental phase shifts in 



4 Tabic |II] indicates that the measured energies E agree better 
with the prediction of the continuum dispersion relation E^^t 



than with the prediction Ef'l' based on free lattice theory. This 
may be due to the smearing of link and quark operators which 
improves rotational invariancc properties. 
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8-13 


1.31/4 
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7-16 


10.1/8 


am K = 0.34660(86) 






(0,0,1) 
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9-16 


4.22/6 


0.5236(11) 


0.52376(58) 


0.51724(58) 


(1,1,0) 


4 




8-14 


0.89/5 


0.6516(24) 


0.65463(46) 


0.64148(45) 



TABLE II: The ground state pion and kaon energies extracted for three momenta: E is extracted from the vari ational a nalysis 
using the chosen interpolator sets. They are compared to analytic expectations in the continuum E^ont — \f m2 + P 2 an d i n 
the free lattice theory Ef a l'a = cosh -1 [cosh(ma) + 2 £\ sin 2 (ip;a)] where the error comes solely from arrin,K above. 



Fig. [2] The phase shifts S^ 2 ' 3 ^ 2 for ^fs ~ + ran near 
threshold are omitted from Fig. [5] and are expressed in 
terms of the scattering length below. The lattice values 
of the phase shifts presented 5 in Fig. [5] apply to y/s quite 
far away from the threshold and we expect that they are 
not significantly influenced by the exact position of the 
threshold, which is at ~ 140 + 500 MeV in experiment 
and at ~ 266 + 552 MeV in our lattice simulation. 

The dependence of the energy levels on the interpolator 
choice is summarized in Fig. [3J We can clearly identify 
the relation between particular levels and the correspond- 
ing meson-meson interpolators, which will be detailed for 
each channel below. 



1. Kn in s wave, 1 — 1/2 

The experimental phase in Fig.[2]is positive and rather 
slowly rising in this attractive channel. The scalar reso- 
nance Xg(800) or k is controversial as it does not render 
a typical Breit-Wigner shape with S = 90° at the posi- 
tion of the resonance. The experimental phase shift does 
not reach S ~ 90° before sfs ~ 1.3 GeV which is in the 
vicinity of the i^(1430). 

In order to understand our lattice spectrum in this 
channel, we first compare the effective energies resulting 
from different subsets of interpolators 0\ q 4 and 05 /A g 
(|A6[) which are plotted in Fig. [3] The horizontal broken 
lines show the energies of the noninteracting scattering 
states K(n)n(— n). The dotted-dashed line corresponds 
to the noninteracting scattering state K*(0)p(0) with the 
energy Ea = m p a + mK*a ~ 0.51 + 0.57= 1.08 (to_r-* can 
be read off from Table IIII1 while m p is provided in Ref. 
[H). The energy of the ground and the first excited states 
is robust to the choice of the interpolator basis as long 
as meson-meson interpolators are included. Using just 
qq interpolators 0\ q 4 gives a much noisier ground state 
and higher energy for the first excited state (right plot 



in Fig. [3]). If we omit Oq from the basis (|A6|) . the state 
K(l)ir(— 1) disappears from the spectrum in Fig. [3] as 
expected. So the relatively noisy level n — 3 (green circles 
at Ea ~ 0.95) can be identified with the back-to-back 
momentum scattering state K(1)tt(— 1). If we omit Oj ~ 
K*(0)p(0), the corresponding state disappears from the 
spectrum, so level n = 4 (blue left-facing triangles) can be 
identified with K *(0)p(0). Removal of O7 has little effect 
on the nearby levels (the resulting -©1,2.3 agree whether 
O7 is the basis or not): thus this interpolator appears to 
be weakly coupled to the system. 

Before presenting the values of the resulting phase 
shifts, we need to discuss up to which E = yfs the Kir 
scattering is elastic on our lattice, since relation © rig- 
orously applies only for elastic scattering. Kir in s wave 
has I{J P ) = 7}(0 + ), which allows the low-lying scatter- 
ing states K*(0)p(0), K(0)r](0) and K(0)r/(0). Since 
our simulation has only dynamical u/d quarks, only the 
heavier 6 772 can occur as an intermediate state in our 
simulation, and we are safe from inelasticity related to 
if (0)77(0). The state if (0)772(0) is expected to be heavier 
and we do not incorporate the corresponding interpola- 
tor, so we expect that our levels n < 4 are not affected by 
that. Our level n = 4 corresponds to the state K *(0)yo(0), 
so the inelastic threshold definitely opens at E±a ~ 1.08, 
while our levels n < 4 are believed to be in the elastic 
regime. 

The effective energies and the final physics results for 
the three levels n < 3 in the elastic regime are presented 
in Fig. [1] and Table IIII1 while the extracted phases are 
compared to experiment in Fig. [2] The phase shift for 
the lowest level near threshold is imaginary and will be 
presented in terms of the scattering length below. The 
phase shifts in Fig. [2] arise from n — 2,3 and agree well 
with the experimental phase shifts. Note that both phase 
shifts are for yfs quite far away from threshold, so the 
exact position of the threshold m n + mx in the lattice 
simulation is not expected to be of major importance 



J Except for lowest level in p wave 1 = 1/2 scattering that is due 
to K* and will be discussed in detail below. 



6 The one that is not the Goldstone boson in SU(2), but is lifted 
by the chiral anomaly. 
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FIG. 1: Effective energies E{i)a of the lowest eigenvalues for the interpolator choices listed in Table Hill together with the 
resulting energies obtained with 1-exponential or 2-exponential fits. The horizontal broken lines show the energies E = Ex + E n 
of the noninteracting scattering states K (n)ir(— n) as measured on our lattice; K{n)n(—n) corresponds to the scattering state 
with p* = y/n^. Note that there is no A(0)7r(0) scattering state for p wave. Black and green circles correspond to the shifted 
scattering states, while the red stars and pink crosses correspond to additional states related with resonances. 



here. 

The level n = 2 corresponds to a phase shift value 
close to 90°, and we attribute this to the vicinity of the 
A'q(1430) scalar resonance. Since 6 for n — 3 is already 
equal to 180° within error bars, we cannot provide a re- 
liable estimate for the resonance mass and width of the 

1/2 

Kq(1A30). Ignoring the huge error bar on p ' (n = 3), 

the central values of p~Q 2 (n — 2,3) lead via a linear 
fit (fT3"|) to the resonance mass ~ 1.25 GeV and 
ryiat ^ q g7 Q e v 2 . This rough estimate of m 1 ^ is within 
T ex P = 270(±80) MeV from m^ p = 1425(±50) McV, 
while -f lat is also not far away from j exp ~ 0.83 GeV 2 
derived from T ex P(Kl 430 -> Kir) = (0.93 ± 0.10) • T exp . 



A reliable lattice determination of these two parameters 
in a future simulation will require an accurate value of S 
near K(l)ir(—1) which is not 180° within error bar. 

The broad resonance AT* (800) (or k) does not lead to 
an additional energy state in our lattice spectrum. This 
is in agreement with expectation since the experimental 
phase shift never reaches S ~ 90° near i/s ~ m K . Similar 
to the 7T7T s wave various analyses assume, that there 
is a broad A'* (800) (or k) resonance hidden behind the 
slowly rising phase shift below 1 GeV; it is associated 
with a resonance pole (in the second sheet) quite distant 
from the real axis. In order to understand that we cannot 
expect an additional level due to k; we plot the analytic 
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8 




n 






range 


type 




[GeV] 








[degrees] 


1=1/2 


1 


01,4,5,6,7,8 


4 


11-16 


1 exp. 


0.4768(28) 


0.7593(45) 


5.8/4 


i 0.0889(31) 


0.409(58) 


i 28.2(5.8) 


s wave 


2 


01, 4,5,6, 7,8 


2 


3-13 


2 exp. 


0.777(13) 


1.237(21) 


9.6/7 


0.2835(87) 


0.051(35) 


82.0(5.3) 




3 


01,2,5,6,8 


4 


6-11 


1 exp. 


0.980(45) 


1.561(72) 


0.42/4 


0.410(26) 


-1.3(3.2) [*] 


162(28) 


1=3/2 


1 


05,6,7,8 


4 


10-16 


1 exp 


0.5323(29) 


0.8478(46) 


0.018/5 f 


0.0653(52) 


-1.67(22) 


-4.21(89) 


s wave 


2 


05,6,8 


4 


8-13 


1 exp. 


0.9979(74) 


1.589(12) 


9.1/4 


0.4208(43) 


-0.76(15) 


-29.1(4.9) 


1=1/2 


1 


01,2,3,6 


4 


8-16 


1 exp. 


0.5749(19) 


0.9156(30) 


10.7/7 


0.1225(21) 


-0.0091(1) 


160.61(73) 


p wave 


2 


01,2,3,6 


4 


8-12 


1 exp. 


0.9558(44) 


1.5223(70) 


0.83/3 


0.3958(26) 


-1.2(1.0) 


177.0(2.6) 




3 


01,2,3,5,6 


4 


6-10 


1 exp. 


1.080(11) 


1.720(17) 


0.32/3 


0.4686(65) 


-0.026(60) 


93.5(7.9) 




4 


01,2,3,5,6 


4 


6-10 


1 exp. 


1.141(18) 


1.817(28) 


1.2/3 


0.503(10) 


0.33(14) 


53(11) 


1=3/2 


1 


06 


/ 


8-14 


1 exp. 


0.9653(31) 


1.5356(48) 


4.5/5 


0.4015(18) 


-0.443(91) 


-8.6(1.8) 


p wave 

























TABLE III: Final results in four channels of Kir scattering: s wave and p wave with 1 = 1/2 and I = 3/2. The total momentum 
of the Kir system is + px = in our simulation. For each Kir eigenstate, we present the energy E = y/s, the momentum 
p* = |Ptt| = |pk| (0! the resulting scattering phase shift S\ ((5} and the quantity p\ defined in Eq. ([§]). The "interpol" column 
indicates which interpolators for the s wave [see Eqs. (|A6I IA7[) ] and for the p wave [see Eqs. (|A8I IA9|> ] are taken as our final 
choice in the variational basis. All fits are correlated with given x 2 (with exception of level n = 1 for 7 = 3/2 in s wave, marked 
by t). Note that the value of pi (HJ has a huge error bar when S is 0° or 180° within error bar, as marked by [*]. Phase shifts 
are determined up to multiples of 180° from Eq. JSJ. 



8 = atan(<77r 3 / 2 /Zoo) [see Eq. @] and the experimental 
phase 5^ 2 as a function of q 2 = (p*L/2tt) 2 with our 
L = 16a in Fig. 2] The energy states on the lattice are 
expected at those q 2 where analytic and experimental 
phases cross if m Vj K is physical. Since our tti^ k are 
not physical, the crossings and energy levels are slightly 
shifted in the actual simulation. The two crossings in 
Fig. 0] correspond to the levels n = 2, 3 in Tabic UTTl while 
there is an additional crossing for n = 1 below threshold 
at imaginary 5 which is not plotted. Obviously, there 
is no crossing near yfs ~ m K , since the lowest crossing 
in the plot 7 at q 2 ~ 0.7 corresponds to y/s well above 1 
GeV. This indicates that no energy level is to be expected 
near yfs ~ m K , just like confirmed by our simulation. We 
emphasize that the absence of an additional energy level 
near yfs ~ m K in our simulation does not contradict the 
possible presence of the k pole in the second sheet; note 
that several analytical studies (for example Ref. [H], Hfl) 
recover the experimental phase shift and do find the pole. 
The lattice simulation is restricted to real s and does not 
have direct access to look for poles as a function complex 
s (just like experiment). So our conclusion is that we 
qualitatively agree with the experimental phase shift in 
this channel, but we cannot conclude whether the k pole 
exists or not. 

Finally, we turn our attention to the lowest level 
slightly below threshold, which has small and imaginary 
p* and 6. It allows the extraction of the scattering length 



from p^ 2 (n = 1) in Table UTTl according to (fT0|) 8 

a /2 = 5.13 ± 0.73 a = 0.636 ± 0.090 fin (16) 

1/2 

^2— = 17.9 ± 2.5 GeV~ 2 at ~ 266 MeV 

where m^- i7r from the simulation were inserted to the re- 
duced mass [J, Kir- In Fig. [5] we compare the values of the 

I—\ /2 

ratio a / hkit with other dynamical lattice simula- 
tions. We choose to present this ratio as it is independent 
of m K ^ in LOChPT 

1/2 3/2 -. 

^ = -2 ^ = -±r 2 [1 + O(ml)} , (17) 

with .F^ ~ 0.13 GeV. This leading-order prediction is 
also shown in Fig. [S] together with the result from ChPT 
at 0(p 4 ) and the result using a Roy-Steiner analysis 
phi . We are not able to perform the extrapolation of our 

a o _1 ^ 2 t° physical pion mass as we calculated its value 
at only one m^. 

As indicated above, we find only one energy level 
below E = 1 GeV in our present simulation: this 
level is related to K(Q)tt(Q), while k does not lead to 
an additional energy level. In our previous simula- 
tion (87l . HH , we employed only four-quark interpolators 



8 Only the first term on the right-hand side of Eq. dlOI I is used to 
7 The crossing due to level n = 2 takes place at q 2 = 0.52(3) in determine Oq 2 as the second term is much smaller due to small 

the actual simulation. p* , and can safely be neglected. 
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FIG. 2: The extracted K ty scattering phase shifts 5j in all four channels / = 0, 1 and / = 1/2, 3/2. The phase shifts are 
shown as a function of the Kn invariant mass i/s = Mktt = y (Pn + Pk) 2 - Our results (red circles) apply for m w ~ 266 MeV 
and rriK — 552 MeV in our lattice simulation. In addition to the phases provided in four plots, we also extract the values of 
ji/ 2 . s/ 2 near threshold ^/s — m-n +iriK, but these are provided in the form of the scattering length in the main text (as they are 
particularly sensitive to m^x). Our lattice results are compared to the experimental elastic phase shifts from Estabrooks (black 
pluses) [llj] and Aston (blue stars) [l^. Dark green crosses represent measured phase shifts by Aston [T3| which correspond to 
an almost elastic amplitude T/, i.e., 0.85 < |2T/ — i\ < 1.15 (see Sec. II A|) . Lattice phase shifts are determined up to multiples 
of 180° from Eq. ©. 



^2 x q(x)q(x)q(x)q(x) with 5 different color and Dirac 
structures, which all had 1=1/2 and J p = + . The 
necessary contractions are shown in Fig. [7^, and [7J:. We 
omitted 9 the backtracking box contraction [7J; in Refs. 
[87l Hl| and the simulation rendered an additional state 
near i^(0)7r(0) which was attributed to k with a sizable 
tetraquark Fock component in Refs. [13, 113 • The effect 
of the box contraction [7J: in our present simulation is 



The first reason for this omission was the numerical cost. The 
second reason had physical motivation of artificially prohibiting 
the mixing qqqq — ¥ qq — ¥ qqqq, so that a qqqq Fock component 
could be attributed to the resulting state. 



shown in Fig. [6l where the spectrum is calculated using 
only (qq)(qq) interpolators ©5,7,8 [see E q (IA6[) ] ; these 
are similar to ©1.2,3 used in Refs. [87l l88|. There is 
only one energy state below E = 1 GeV when all nec- 
essary contractions ([7Jt and [7J:) are incorporated, which 
agrees with the result in Fig. [T]and with our conclusions 
above. However, if the backtracking box contraction^: is 
neglected, an additional energy level near i^(0)7r(0) ap- 
pears. A proper quantum field theory treatment requires 
incorporation of all Wick contractions, so the additional 
level seems to be an artifact of the approximation used 
in [13, [13 ■ This interesting observation may be fruitful 
in trying to understand the physics of light scalar mesons 
in future explorations. 
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FIG. 3: Effective energies E(t) a of the lowest eigenvalues for different interpolator choices in the correlation matrix. The 
horizontal broken lines show the energies E = Ek + E^ of the noninteracting scattering states K(n)n(—n) as measured on our 
lattice with p* = 2zl ^Jn. Note that there is no K(0)ir(0) and K*(0)p(0) scattering state for the p wave. Red stars and pink 
crosses correspond to states related to resonances; other levels are related to scattering states, (a) s wave, 1=1/2: the first four 
choices incorporate also meson-meson interpolators Ol^.g; the fifth choice incorporates just qq interpolators 0\ q 4 [see Eq. 
(|A6|) ]. (b) s wave, 1=3/2: various choices of meson-meson interpolators Ol^g [see Eq. (|A7[) ]. ( c ) P wave, 1=1/2: levels with 
and without meson-meson interpolator Oa [see Eq. (|A8|) ] in the basis. 



2. Kir in s wave, I = 3/2 

This K + ir + channel is repulsive, and the experimental 
phase shift in Fig. [2] is negative and slowly rising. A 
negative phase shift in Table IIIII is also observed for the 
lowest two states in our lattice simulation, since they are 
clearly above noninteracting K(0)ir(0) and K (!)%(— 1) in 
Fig. [0 The value of 5 ~ -30° at y/s ~ 1.6 GcV agrees 
nicely the experiment. 



We do not see any additional state between these two 
levels, which agrees with the fact that resonances have 
not been experimentally observed in this repulsive chan- 
nel. 

Investigation of the spectrum with different interpola- 
tor choices in Fig. [3] indicates that Oq [see Eq. (|A7[) ] is 
responsible for level n = 2, and O-j ~ K*(0)p(0) is re- 
sponsible for the level n — 3. The Krj and K rj do not 
contribute to the 1 = 3/2 channel due to isospin, so the 
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4: The analytic phase 5 = atan[7r 3//2 g/Zoo(g 2 )] (see 
£jq. @ ) as a function of q 2 . The experimental phase 5^^ 2 
[ill [l2j| is also given as a function of q 2 = (p*L/2-rr) 2 , where 
L — 16a is our lattice size and p* is the center-of-momentum 
frame of n or K in the experiment. The energy levels in a 
lattice simulation arise for the values of q 2 where the analytic 
and the "experimental" phases cross. Since the simulation 
is not performed at physical itik^ the crossings are slightly 
shifted in the actual simulation. 



lattice data and the experimental data are completely 
clastic up to rather high yfs = niK* + m p . 

The lowest level slightly above threshold allows the 
extraction of the scattering length from p^ 2 (n = 1) in 
Table Hill according to (fTO)) 



7=3/2 
1=3/2 



= -1.13 ± 0.15 a = -0.140 ± 0.018 fm 



(18) 



-3.94 ± 0.52 GcV^ at ~ 266 McV . 



This scattering length is compared with other lattice and 
continuum results in Fig. [5J 

Due to the smooth behavior of the phase shift observed 
in experiment [ll| , we attempt to estimate also the effec- 

3/2 

tive range r by employing the effective range formula 

(fTU)) . From the values of p^ 2 , p* and -y/s for the levels 
n = 1, 2 in Table IIII1 we extract 



a 3 Q /2 = -1.12 ± 0.15 a = -0.139 ± 0.018 fm (19) 
r Z=3/2 = 15± 2.0 a = 0.19 ±0.25 fm 

at our m n , which indicates that the dependence of p* cot S 
(|10p on p* is small (zero within errors) up to p* ~ 
0.67 GeV. Our I = 3/2 phase shift is therefore dictated 
by the scattering length for p* as high as p* ~ 0.67 GeV. 
The experimental effective range at physical pion mass is 
r 3 Q /2 = -0.346(±0.060) fm [ll|. 



3. Kn in p wave, 1—1/2 

The elastic region is dominated by a vanilla-style res- 
onance: the A* (892) with a width of w 50 MeV from 
experiments. Experiments indicate further resonances 
if* (1410) and A'* (1680), where the first one is not es- 
tablished in all experiments. 

The spectrum on our lattice is shown in Fig. [TJ and 
all these levels are expected to be in the elastic regime, 
as discussed below. The scattering level A(l)7r(— 1) is 
seen only if the meson-meson interpolator Oq [see Eq. 
(|A8[) ] is taken in the interpolator basis, while qq interpo- 
lators alone do not render it (see Fig. [3]). The other 
three levels (red stars and pink crosses) appear away 
from the noninteracting scattering states K(n)ir(—n) 
and are candidates to be related with the resonances 
A*(892), IT* (1410), A*(1680). 

The ground state is due to the A"* (892) and gives a 
rather high S ~ 160° at y/s ~ 915 MeV ~ m K * (see 
Fig. [5]). This is not surprising since the phase should be 
rising very steeply as expected from the narrow width 
r' at ~ {p*iatlP*ex P ? Texp - (0.19/0. 29) 3 • 50 MeV ~ 
14 MeV derived from Y exp and assuming the same cou- 
pling 7 [see Eq. ([T2"j) ] in both cases. Since we have only 
one value of the phase shift near \/s ~ mx * , we can- 
not rigorously extract the A* (892) mass or width, but 
we expect that the phase shift would pass 90° at about 



lat 



+ yfiT lat cot S) 1 ' 2 ~ 896 MeV (O, where the 
derived width T' at quoted above was assumed. Extrac- 
tion of the A* (892) width in future lattice simulations 
will be possible only if two phase shifts are extracted in 
close vicinity of \/s ~ ma* ■ Simulations at nonzero-total 
momentum and the relevant phase shift formulae for the 
p wave [12, [7l,[73 might come to the rescue in this case. 

The third and fourth level at E 3 = 1.720(17) GeV and 
E4 = 1.817(28) GeV arc most probably related to the 
wide resonances A* (1410) , A* (1680) which are the only 
p wave resonances which appear between 1 and 2 GeV in 
experiment. However, the level A3 appears too high in 
comparison to the A* (1410) even if one takes into consid- 
eration its large experimental width Y exp = 232(21) MeV 
and our unphysical rrin.K- Notice, however, that in a 
small box and with unphysically heavy u/d quarks the 
situation is quite different from experiment, where both 
the A* (1410) and the A* (1680) have a sizable branching 
ratio into the A*(892)7r and pK channels. The inelastic 
threshold in this J p = 1~ channel opens at K*(l)ir(— 1) 
or A(l)p(— 1) in p wave. This is at rather high A ~ 1.9 
GeV, so all levels (except possibly n = 4) are in the elas- 
tic regime in our simulation. Therefore, our situation is 
somewhat unphysical in this case, as we can only consider 
elastic scattering in the A-7T channel. Our results are con- 
sistent with the observations in the simulation [89| , where 
the energy level associated with the A* (1410) is also ob- 
served substantially higher than the physical state. No- 
tice that for our kinematics, one expects a phase shift 
which is monotonically increasing in the vicinity of the 
levels A3 and A4. This is not in conflict with Fig. [21 where 
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FIG. 5: The s wave scattering lengths clq for / = 1/2 and 3/2 expressed as clq/ ' Hk-k which does not depend on at LOChPT. 
The result from the present simulation at a sing le ~ 266 MeV is compared with results from other dynamical simulations 
[HQ and with LOChPT, ChPT at 0(p 4 ) [53 and a Roy-Steiner approach Q. 



all contractions 



ignoring box diagram 



1.2 



3 

^0.8 



0.6 



0.4 



10 

t 



15 



5<j>c|> 



W © © © © © © ©©© © © © ©■ 



10 

r 



15 



FIG. 6: Effective energies in the s wave J = 1/2 channel 
obtained from a 3 x 3 correlation matrix using only meson- 
meson interpolators Os,7,8 (see l|A6|l ) and to — 4. Left: using 
all necessary contractions C direct + ^C crossed - \C hox listed in 
Appendix B and illustrated in Figs. [7K and[7h. Right: using 
only C direct + lC crossed (Fig. Hi) and omitting the backtrack- 
ing box contraction in Fig. [TE 



the phase is restricted to < S < 180°, since the phase 
obtained from Eq. ([5]) is undetermined up to multiples 
of 180°. 



On the lattice, we extract the phase from a single 



level given by the correlator (0 6 (£)|0g(O)) 



-El 



see 



4- K-k in p wave, I — 3/2 



Eq. (|A9|) ]. It appears slightly above the noninteract- 
ing K(l)ir(— 1) and renders a negative phase S lat (y/s ~ 
1.5 GeV) = -8.6° ± 1.8°. This is in qualitative agree- 
ment with the experimental phase which is negative and 
does not exceed —10° for ^/s as large as 1.8 GeV (see 
Fig. [2]), although both phases do not agree within errors. 



IV. SUMMARY AND OUTLOOK 

We simulated Kit scattering with lattice QCD and ex- 
tracted the elastic phase shifts b\ in s wave and p wave 
for / = 1/2,3/2 at several values of the Kir invariant 
mass v^s. We used a single lattice QCD ensemble of size 
L ~ 2 fm with dynamical u and d quarks. Our results for 
phase shifts and scattering lengths apply for the values 
of m T ~ 266 MeV and m K ~ 552 MeV. The total three- 
momentum of the Kit system is zero in our simulation. 

First, we extracted the energy levels of the Kir sys- 
tem on our finite lattice. In all channels, we observe 
the expected K{n)n(— n) scattering state levels, which 
are shifted relative to the nonintcractive case due to 
the interaction. In both attractive I = 1/2 channels, 
we observe additional levels which are related to reso- 
nances A^(1430) in s wave and A'* (892), K* (1400) and 
iT*(1680) in p wave. 

The phase shifts are extracted from the energy levels 
using Luscher's method. They are compared to the ex- 
perimental phase shifts in Fig. [2] and exhibit qualitative 
agreement in all four channels: 



The experimental phase in this repulsive channel is 
negative and very small. 



• s wave, 1=1/2: The phase is positive and yields 
the scattering length Oq -1 ^ 2 = 0.636(90) fm at our 
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m n . Our first excited state is observed with <5 ~ 90° 
near Aq(1430), which implies that the phase does 
not reach 90° below s/s = 1 GeV in our simula- 
tion. This agrees with the experimental finding 
that the phase is not 90° near a/s ~ m K and that 
the controversial k resonance cannot be described 
by a conventional Brcit-Wigncr shape. 

• s wave, 7 = 3/2: The phase is negative and reaches 
(5q _3//2 ~ —30° at y/s ~ 1.6 GeV in agreement with 

experiment. We extract a^~ 3 ^ 2 = —0.140(18) fm 
and we find very mild dependence of y/s cot d on p* 
up to p* = 0.67 GeV. 

• p wave, 1=1/2: Our spectrum and the phases 
favor the existence of three resonances if* (892), 
if* (1410) and A"* (1680) below 2 GeV, but our en- 
ergy level for A"* (1410) is higher than the experi- 
mental one. We did not extract the A* (892) width 
as we have only one energy level in the vicinity of 
this well-established narrow resonance. 

• p wave, 1 = 3/2: The phase is negative and very 
small, which is observed also in experiment. 

The extraction of the widths for the resonances in the 
i = 1/2 channels is left for future simulation. This 
will require at least two values of the phases S(y/s) in 
the vicinity y/s ~ mji ± Tji of each resonance, which is 
particularly challenging for a narrow resonance like the 
A* (892) . This might be possible to achieve for the p wave 
using simulations with P = + ^ which would 
provide additional values of phases 5g = i at s = E 2 — P 2 . 
The relations which allow an extraction of <5i=i from the 
energies in this case are derived in Refs. [H, IzE l77| . 
while sample interpolators arc explicitly listed in Ref. 
[76j. This will be much more challenging for the s wave 
since Se=o is always mixed with St=i in the phase shift re- 
lations for the A\ irreducible representation when P^0 
and ^ rriK [H, [7514771 . The extraction of the s-wave 
phase shifts is therefore more reliable with the present 
simulation at P = and calls for similar simulations at 
different lattice sizes. 
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Appendix A: Interpolators 

1. Single pion and kaon interpolators 

For the single pion or single kaon sectors we have 6 
interpolators, using three smearing widths [see Eq. (| 15[) ] 
for each of the two Dirac structures, 

0? ype , s (p,t) = ^3,(i)r fwt e'Pu,(i) ) 

X 

0? ype , s (p,t) = J2Mx)r tyP ee ip *u s (x) , (Al) 

x 

Ti = 75, r 2 = 757t, s = n,m,w. 

Here we determine the energy levels for different values 
of the total momentum p in order to study the dispersion 
relation. 



2. Interpolators for the Kn system with P = 

For the Kit system we use qq and meson-meson (Ktt, 
pK*, Kicli) interpolators with appropriate quantum 
numbers. The interpolators for all four channels are 
listed in the following subsections. 

The meson-meson interpolators arc expressed as prod- 
ucts of two meson currents, where each meson current 
M (p) has momentum p projected to or ^e, 

Ai(p) = ^g 1 (x,t)re J P x (?2 (x,t) 



We use the following flavor combinations and T matrices 



for M(p) 








7T+ 


= dj 5 u , 


7T° = 


i/f (M75U - dj 5 d) 


A+ 


= S75M , 


A° = 


sj 5 d 


(p + )i 


= djiU , 


(A = 


= y/% (HjiU - djid) 




= SJiU , 


(A* ), 


= sjtd 




= drfn&u , 




= (u7i75M - dji^d) 




= S7i7 5 u , 


(*?)< 


= sji"f 5 d 



where i = x,y,z refers to the three spatial directions. 
Each quark q s is smeared according to Eq. (|15[) with the 
number of eigenvectors N v , provided for the interpolators 
below. 

We combine these meson currents into the I = 1/2 
combination 

\I,h) = \hi) = Vl^ + V / |A' 7r+ (A3) 
or the 1 = 3/2 combination 

|/,/ 3 > = U, i)=K+ir+ (A4) 
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and analogously for (p,K*) and (ax,Kx) pairs which 
carry different J p and the same isospin. 

Our quark-antiquark interpolators below contain also 
a covariant derivative, defined as 

^i(x,y) =C/ i (x,0)<5 x+i>y -C//(x-i,0)J x _ i)y . (A5) 

It acts on the spatial and color indices and leaves time 
and Dirac indices intact. 



3. Kir in s wave, 1—1/2 

For the k channel we employ 4 quark-antiquark inter- 
polators and 4 meson-meson interpolators in the varia- 
tional basis. Interpolators Oi_5 are built using q m with 
N v = 64 eigenvectors, while interpolators ©6-8 take q w 
with 32 eigenvectors due to the sizable numerical cost re- 
lated to them. The quark-antiquark interpolators 0i_4 
differ in Dirac and color structure. The interpolator O5 is 
a Kit interpolator where both pseudoscalars are at rest, 
whereas for Oq they have oppositely oriented unit mo- 
mentum, summed over all spatial directions. Finally Oq 
and Or are pK* and a\K\ at rest. These two are in the 
inelastic region and their (ir)relevancc is discussed in Sec. 
IIIIB II So we compute an 8 x 8 correlation matrix with 



O x = ^2 s(x)u(x) , 

X 

2 = ^2 s( x ) V * u(x) , 

x.v 

x,i 

5 = J^K + (0)tt°(0) + JzK°(0)tt + (0) 



(A6) 



06 = £ [ypr+( Pi )7r (- Pi ) + vTtf°(Pi)7r + (-Pi)] 

i 

+ (Pi <-> -Pi) , 

i 

Os = E [^,+(0) «?,(()) + y|K 1 °(0) a + t (0)] , 



and the sum on i runs over i = x,y, z. The momenta P i 
are given in Eq. (|A2[) . 



4. if7r in s wave, 1 — 3/2 

For the exotic I — | c hannel we use the corresponding 
interpolators O5-8 (|A6[I and the same choice of smear- 



ings, just a different isospin projection 

O 5 =/v + (0)7r+(0), (A7) 
O e = K + ( Pi )7r + (- Pi ) + K + (- Pi )7r + ( Pi ) , 

i 

(9 7 = ]Tif*+(0)p+(0), 

i 

O 8 =^K 1 +(0)a+ i (0). 

i 

The naming scheme is kept analogous to Eq. (|A6I) . 

5. Kit in p wave, I = 1/2 

For the X* channel we employ quark-antiquark inter- 
polators 0i_5 and one kaon-pion interpolator Oq, which 
are all built using q n with N v = 96 



(A8) 



2 ,r 



^2 s(x) 7*7, u(x) , 

X 

^2 s(x) V j 7i Vj u(a 



u(x) , 



m(x) , 



Oh,i = ^2 Ci i k s ( x hjl5 \ 
- (Pi ^> -Pi) • 



O 



Here, the open polarization index is i = x, y, z, and 
we average the resulting correlation matrices over three 
polarizations. The linear combinations of the deriva- 
tives in 04,5 render good C parity in the SU(3) flavor 
limit. Although C is not a good quantum number due 
to m s ^ m Ut d, such a combination is advantageous as 
discussed, e.g., in Ref. (89j . 



6. Kit in p wave, / = 3/2 

For the I = 3/2 p wave channel, we use only Oq from 
Eq. (|A8[) with appropriate choice of flavor 

O e ,i = A+( Pl )7r + (- Pl ) - K+(- Pi )Tr + { Pi ) , (A9) 

where i = x,y,z is the polarization index. As there is 
just one interpolator in this case, we present its complete 
expression for convenience 

6>i = J2 [e i(PtXl " PiX2) ^(i,x 1 ) 75W „(i,x 1 ) 

X1,X 2 

X d„(i,X 2 )75Wn(*,X2)] - ( P j <-> -Pi) ■ 

Other employed meson-meson interpolators can be ex- 
pressed in terms of the quark fields in an analogous way. 
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FIG. 7: Contractions for our correlators with qq and meson- 
meson interpolators. Only (a) appears for / = 3/2, while all 
these contractions appear for / = 1/2. The "backtracking" 
contractions (c) and (d) require an all-to-all method. 



Appendix B: Wick contractions 

Here, we provide the expressions for the 1=1/2 and 
7 = 3/2 contractions. In the distillation method, they 
are expressed in terms of the perambulators for the light 
quark t u and the strange quark t s , as well as <f> matrices 
which depend on the shape functions J 7 . We use exactly 
the same definitions of these quantities as in Sec. IIC 
and Appendix A of Ref. [l| , so we omit the definitions 
here. 

Any annihilation operator 0(t) above can be expressed 
in terms of 

MM (t) = y§ s(t)T 2 ^{p 2 )u(t) (Bl) 

x ^\ [u(i)ri^(pi)«(t) - d{t)Y x T{p x )d{t)\ 

+ yfl s{t)T 2 F{p 2 )d{t) d{t)T 1 T{p 1 )u{t) , 
O qq (t) = s(t) T T a (P) u(t) , 

and creation operators 0(t') at source can be expressed 
as 

0\ IM (t') = Cum [y/l u(t')T' 2 F(-p 2 )s(t') (B2) 

x yf\ w^h-p'iW) - dWr'^i-tiW)] 

+ yfi d(t')T' 2 F(- P ' 2 w) uitn^-tiw)] , 

Ot q (t') = C qq u{t>) T' T (-P) s(t') . 
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